Enhanced Non-local Power Corrections to the B — ^ Xsj Decay Rate 
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A new class of enhanced non-perturbative corrections to the inclusive B —> Xs^ decay rate is iden- 
tified, which contribute first at order A/mt in the heavy-quark expansion and cannot be described 
using a local operator product expansion. Instead, these effects are described in terms of hadronic 
matrix elements of non-local operators with component fields separated by light-like distances. They 
contribute to the high-energy part of the photon-energy spectrum but do not reduce to local oper- 
ators when an integral over energy is taken to obtain the total inclusive decay rate. The dominant 
corrections depend on the flavor of the B-meson spectator quark and are described by tri-local four- 
quark operators. Their contribution is estimated using the vacuum insertion approximation. The 
corresponding uncertainty in the total decay rate is found to be at the few percent level. This new 
effect accounts for the leading contribution to the rate difference between B~ and B° mesons. 
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I. INTRODUCTION 

Precision studies of inclusive i?-nieson decays are a cor- 
nerstone of quark flavor physics. Detailed measurements 
of various kinematical distributions in the semileptonic 
decays B —^ XI D, when combined with elaborate the- 
oretical calculations, provide the currently most precise 
measurements of the elements \ Vcb\ and \ Vub\ of the quark 
mixing matrix (see [1] for a comprehensive recent anal- 
ysis of charmless inclusive decays). Studies of the rare 
decays B Xg^ and B ^ Xg I'^l' allow sensitive tests 
of the flavor sector and provide constraints on extensions 
of the Standard Model. 

The theoretical description of inclusive _B-meson decay 
rates is based on the operator product expansion (OPE) 
[2,3], by which total decay rates can be expressed in terms 
of forward _B-meson matrix elements of local operators. 
Only two non-trivial matrix elements appear up to or- 
der {K/mbY in the expansion, one of which can be ex- 
tracted from spectroscopy. The OPE breaks down when 
one tries to calculate differential inclusive decay distri- 
butions near phase-space boundaries. A twist expansion 
involving forward matrix elements of non-local light-cone 
operators (so-called shape functions) is then required to 
properly account for non-perturbative effects [4,5]. Re- 
cently, these non-local structures have been analyzed sys- 
tematically beyond the leading order in A /mi, [6-8]. It is 
generally believed that the non-local operators reduce to 
local ones when the differential decay distributions are 
integrated over all of phase-space. Here we show that 
this is not always the case. 

A precise control of hadronic power corrections is par- 
ticularly important in the case of the inclusive radiative 
decay B Xs^, which is the prototype of all flavor- 
changing neutral current processes. A significant ef- 
fort is currently underway to complete the calculation of 
the leading-power (in A/mt,) contribution to the decay 
rate at next-to-next-to-leading order in renormalization- 



group improved perturbation theory. This leaves non- 
perturbative power corrections as the potentially largest 
source of theoretical uncertainty. 

It is well-known that in B Xgj decay the OPE 
faces some limitations, which result from the fact that 
the photon has a partonic substructure. For instance, 
there exists a contribution to the total decay rate involv- 
ing the interference of the & — *■ 57 transition amplitude 
mediated by the electro-magnetic dipole operator Q^^ 
with the charm-penguin amplitude {b — )■ ccs followed by 
cc —f 75) mediated by the current-current operator Qi 
(see [9] for the definition of the operators in the effective 
weak Hamiltonian) . When the charm-quark is treated 
as a heavy quark (rric ~ rrib), this contribution can be 
expanded in local operators [10-13], and it is believed to 
be a good approximation to keep only the first term in 
this expansion. Its contribution to the total B Xgj 
decay rate can be written as 
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is the leading-order contribution to the decay rate from 
the electro-magnetic dipole operator. The ratio AP/Pyt 
therefore provides an estimate of the relative magnitude 
of the non-perturbative effect. We stress that when the 
scaling ~ Arub is adopted instead of rric ^ rrib, then 
the charm-loop contribution must be described by the 
matrix element of a non-local operator [11-14]. 

It has been noted in [11] that the OPE for B Xg'y 
decay breaks down when one includes diagrams from op- 
erators other than Q7~f, in which the photon couples to 
light quarks. An example of such a contribution, result- 
ing from the decay b sg mediated by the chromo- 
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magnetic dipole operator Qsg followed by photon emis- 
sion from the light partons, was studied in [15]. It was 
argued that the corresponding effect can be estimated in 
terms of the parton fragmentation functions of a quark 
or gluon into a photon. Since this contribution is nu- 
merically very small, it has not received much further 
attention in the literature. A more careful analysis re- 
veals that the correct interpretation of this effect is in 
terms of a subleading shape function [14]. 

In this Letter, we identify and analyze a novel class of 
non-local power corrections to the B —^ X^j decay rate, 
which were not considered before. We argue that they 
can affect the total decay rate at the few percent level, 
and that they give the dominant flavor-specific contribu- 
tion to the rate difference between charged and neutral B 
mesons. The presence of this effect leads to a dominant 
and irreducible source of theoretical uncertainty in the 
prediction for the total B — > Xs^ branching fraction. 

II. NON-LOCAL POWER CORRECTIONS 

Power corrections to the high-energy part of the B 
Xsj photon spectrum can be systematically parameter- 
ized in terms of subleading shape functions defined in 
terms of forward i?-meson matrix elements of non-local 
light-cone string operators [6-8] . Some of these operators 
- the ones considered so far in the literature - reduce to 
local operators when one considers the total decay rate 
(i.e., the integral over the photon spectrum); however, a 
detailed analysis shows that several of them do not [14]. 
Representative diagrams giving rise to such operators are 
depicted in Figure 1. The graphs show different contri- 
butions to the discontinuity of the hadronic tensor W^"^ , 
which determines the B —^ Xs"f photon-energy spectrum 
via the optical theorem. The total decay rate is obtained 
by an integration over the photon energy. In this Letter 
we focus on the two graphs shown in the first row (and 
two mirror graphs, in which the order of the weak ver- 
tices is interchanged). They describe the interference of 
the b sj transition amplitude mediated by the electro- 
magnetic dipole operator with the b sg amplitude 
mediated by the chromo-magnetic dipole operator Qgg 
followed by the fragmentation of the gluon into an en- 
ergetic photon and a soft quark-antiquark pair. While 
other diagrams, such as the first graph in the second row 
in the figure, give rise to four-quark operators containing 
strange quarks, the graphs in the first row produce all 
light-quark flavors. We expect that the resulting four- 
quark operators will have a larger overlap with the B- 
meson states and thus give rise to the dominant power 
corrections. For simplicity, we also do not consider loop- 
suppressed effects such as the second graph in the second 
row of the figure. This diagram would match onto a non- 
local operator containing a soft gluon field, which mixes 
with the operators we consider under renormalization. 

The top two diagrams in Figure 1 affect the B ^ Xgj 




FIG. 1. Diagrams representing enhanced non-local power 
corrections to the B — > X^^ photon spectrum, which do not 
reduce to local power corrections to the total decay rate. The 
double lines represent heavy-quark fields hv The vertical 
dashed lines indicate cuts of the relevant propagators. 



photon spectrum in the high-energy region, Ej w mi,/2, 
where it is most accessible to experiment. This is en- 
forced by the fact that the amplitude mediated by the 
insertion of Qgg interferes with the two-body decay am- 
plitude mediated by the insertion of Qj^. The effect can 
therefore not be eliminated using kinematical cuts. We 
find that the corresponding contribution to the total de- 
cay rate can be parameterized in terms of forward B- 
meson matrix elements of tri-local light-cone operators 
(with Cf = 4/3 for = 3 colors): 

C77 Ncrrib J_oo J -00 
X {B\Cf{Oi+02)-{Ti+T2)\B), (1) 

where we have used that the Wilson coefficients Ci are 
real in the Standard Model. The relevant factorization 
scale to use in this result is of order fi^ ^ mi, A. We 
use a mass-independent normalization of meson states, 
such that {B\hyhy\B) = 2. The four-quark operators 
are defined as 

Oi^^Cq hy{0) Tji q{tn) q{sn) T n /i„(0) , 
9 

^2 - II y ^-(0) rfl7±a q{tn) q{sn) jITr K{0) , 
9 

Ti = ^eq hy{0)TRtaq{tn) q{sn) Ffl, ta/i^(0) , 
9 

T2 = ^^hy{0)TRj±c.taq{tn) q{sn)jZTRtahy{0), 
9 

where is the electric charge of the soft light quark in 
units of e, Tr = ^{1 + 75)72 is a right-handed Dirac 
structure, and ta are the generators of color SU(3). Here 
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FIG. 2. Diagram representing an enhanced non-local power 
correction to the total B — > X^j decay rate. 

are the two-component heavy-quark fields defined in 
heavy-quark effective theory, while q arc soft light-quark 
fields (g = u,d,s). The light-quark fields are located 
on the light-cone defined by the direction of the emitted 
photon with momentum = E^n'^ (with v? = 0). The 
non-local operators are made gauge invariant by the in- 
sertion of soft Wilson lines Sn in the fi-direction. The 
Wilson lines are absent in light-cone gauge n ■ A — 0, 
which we adopt implicitly to simplify the notation. 

Evaluating the hadronic matrix elements in (1) us- 
ing a systematic non-perturbative approach is a very 
challenging task. In particular, lattice QCD is unable 
to handle operators with component fields separated by 
light- like distances. Naive dimensional analysis suggests 
that Ar/Fyy ^ {Csg/Crj) nas{A/mb), which could eas- 
ily amount to a 5% correction to the decay rate. In more 
traditional applications of the OPE to inclusive i?-meson 
decays, four-quark operators contribute at order [A/mt,)^ 
in the heavy-quark expansion. The non-local operators 
in (1) lead to enhanced power corrections of order A/mb, 
because the two "vertical" propagators in Figure 1 have 
virtualities of order uibA and so introduce two powers of 
soft scales in the denominator. This mechanism was first 
studied in [16]. The existence of such enhanced power 
corrections in the total decay rate may seem puzzling at 
first sight. Consider the diagram in Figure 2, which rep- 
resents the contribution to the total rate derived from 
the first graph in Figure 1. Before taking the cut indi- 
cated by the vertical dashed line the diagram only re- 
ceives contributions from hard loop momenta p'' ~ mt,, 
and it would thus seem appropriate to shrink all propaga- 
tors to a point. In this case the diagram would contribute 
at order (A/mf,)^ in the heavy-quark expansion, and this 
scaling would appear to be preserved when one takes the 
discontinuity of the diagram, apparently contradicting 
our conclusion. The loophole is that the contribution to 
the total B Xg^ decay rate is not given by the discon- 
tinuity of the loop graph, which would correspond to the 
sum of all three possible cuts, but instead it is given by 
the single cut shown in the figure. 

In order to obtain at least some model estimate of the 
magnitude of the effect in (1) we adopt the vacuum in- 
sertion approximation (VIA), in which the vacuum state 
|0)(0| is inserted between the two light-quark fields inside 
the four-quark operators. This is a crude approximation, 
which however appears to work well in the analysis of 
&-hadron lifetimes [17,18]. The approximation has thus 
been checked for local four-quark operator evaluated be- 



tween i?-meson states. Also, it can be justified using 
large- A'c counting rules. Applications of the VIA to non- 
local operators can be found in [8,19] 

In the present case, the matrix elements of the oper- 
ators O2 and Ti^2 vanish in the VIA, either due to the 
color-octet structure of the quark bilinears (T'1^2) or due 
to the fact that there is no external perpendicular Lorentz 
vector available (O2 and T2). The matrix element of 
Oi can be expressed in terms of the leading light-cone 
distribution amplitude of the _B-meson in position space 
[20,21]. We obtain 

{B\ Or |B)viA = e, ^f{s) [0f (t)]* , 

where now refers to the charge of the light spectator 
quark in the B meson. The integral over the position- 
space distribution amplitude can be evaluated to yield 

-i / ds0f(,s)= / —cj>f{u;) = —, 
J-00 Jo ^ -^B 

where <j>^{uj) is the light-cone distribution amplitude in 
momentum space, and A^^ is the common notation for 
the first inverse moment of this quantity [9]. Numerical 
estimates of As are very uncertain, but typically fall in 
the range between 0.25 and 0.75 GeV [9,20-24]. In the 
VIA, our estimate of the spectator-dependent, non-local 
power corrections then takes the final form 

AFyiA _ CqCsg TTUs A 1_\ /Img 

T77 " Cr^ 2 \ N^J X%mb ' ^ ' 

Recalling that fs ~ ^/^/mB in the heavy-quark limit, we 
indeed recover the scaling behavior anticipated above. 
Numerically, with /i ^ 1.5 GeV as a typical factoriza- 
tion scale and /b ~ 0.215 GeV for the i3-meson decay 
constant (see [25] for a recent determination using un- 
quenched lattice QCD), we obtain 

^._0.8e,(||)^-0.05e.(^)-^ 

where Cq = 2/3 for decays of B~ mesons, while = 
— 1/3 for decays of B^ mesons. For the range of Xb values 
quoted above, the effect is between —2% and —19% times 
Cq. Taking Fyy as an estimate of the total decay rate 
at leading power, this implies that the enhanced power 
corrections to the total, flavor-averaged B Xgj decay 
rate are expected (in the VIA) to be between —0.3% and 
—3%, while these effects induce a flavor-dependent rate 
asymmetry 

T{B- ^ X,7) - F(^° ^ X,7) ^ _ / Ab \ 

T{B^X,-f) ~ %0.5GeVy ' 

(3) 

which could amount to an effect between —2% and —19%. 
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When considering these estimates one should keep in 
mind that the VIA can at best provide a very simple 
model of the effect of the non-local fonr-riuark operators 
in (1). Conservatively, we can therefore not exclude that 
the type of enhanced power corrections identified in this 
Letter could contribute to the total B Xsj decay rate 
at the 5% level. The magnitude of the flavor-specific ef- 
fects studied above could be probed by a measurement of 
the flavor asymmetry (3); but there are other four-quark 
operator contributions with flavor structure bs sb (see e.g. 
the bottom left diagram in Figure 1), whose matrix ele- 
ments vanish in the VIA but could still be significant in 
real QCD. Their contributions are flavor-blind and hence 
not tested by (3). 

The Babar collaboration has measured the flavor- 
dependent rate asymmetry in eq. (3), finding the value 
(1.2 ± 11.6 ± 1.8 ± 4.8)% , where the errors are statisti- 
cal, systematic and due to the production ratio / B~ , 
respectively [26]. The dominant error is statistical and 
therefore likely to decrease when more data is collected. 



III. CONCLUSIONS 

We have identified a new class of enhanced power cor- 
rections to the total inclusive B Xs^ decay rate, which 
cannot be parameterized in terms of matrix elements of 
local operators. These effects are nevertheless "calcula- 
ble" in the sense that they can be expressed in terms of 
subleading shape functions. At tree level, the correspond- 
ing operators are tri-local four-quark operators. While 
local four-quark operators contribute at order {K/mbY 
in the heavy-quark expansion of the total decay rate, 
the effects we have explored arc enhanced by the non- 
local structure of the operators and promoted to the level 
of A/mft corrections. We have identified and estimated 
what we believe are the dominant corrections of this type, 
namely those that match the flavor quantum numbers of 
the external B-meson states. 

Oiiv results imply that a local operator product expan- 
sion for the inclusive B X^j decay rate does not ex- 
ist. Even at first order in A/nib there are hadronic effects 
that can only be accounted for in terms of non-local oper- 
ators. The precise impact of these power corrections will 
be notoriously difficult to estimate using our present com- 
mand of non-perturbative QCD on the light cone. While 
a naive estimate using the vacuum insertion approxima- 
tion suggests that the effects arc at the few percent level, 
we conclude that they are nevertheless a source of sig- 
nificant hadronic uncertainty in the calculation of partial 
or total B — > Xgj decay rates. After the perturbative 
analysis of the decay rate will have been completed, the 
enhanced non-local power corrections will remain as the 
dominant source of theoretical uncertainty. A measure- 
ment of the flavor-dependent asymmetry (3) could help 
to corroborate our numerical estimates of such correc- 
tions. 
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